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Abstract 

We investigate the temporal evolution of a ferromagnetic system of Ising spins evolving under 
Kawasaki dynamics from a random initial condition, in spatial dimensions one and two. We ex- 
amine in detail the asymptotic behaviour of the two-time correlation and response functions. The 
linear response is measured without applying a field, using a recently proposed algorithm. For 
the chain at vanishingly small temperature, we introduce an accelerated dynamics which has the 
virtue of projecting the system into the asymptotic scaling regime. This allows us to revisit criti- 
cally previous works on the behaviour at large time of the two-time autocorrelation and response 
functions. We also analyse the case of the two-dimensional system at criticality. A comparison 
with Glauber dynamics is performed in both dimensionalities, in order to underline the similarities 
and differences in the phenomenology of the two dynamics. 

PACS numbers: 05.70.Ln, 64.60.My, 75.40.Gb 
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I. INTRODUCTION 



The kinetics of ferromagnetic spin systems evolving after an initial quench from a high 
temperature disordered initial condition to a final temperature, equal or below the criti- 
cal temperature, is a well investigated field (see e.g. ref. Q for a review). However, only 
recently has the emphasis been put on the dynamics of two-time quantities, such as the 
correlation and response functions, or the fluctuation-dissipation violation ratio, with the 



aim of quant i: 



tion 



P. 



'ymg^ the distance of the system to equilibrium during its temporal evolu- 



10, 



11 



14L Il5l ll6L 1 1 71. Il8j] (see e.g. [19| for a brief review). 



Most of these studies focus on the non-conserved order parameter case, where at each time 
step a single spin is updated according to the rules of Glauber dynamics. An overall coherent 
picture of this field has by now emerged, though some controversies remain (20, 21, 22I l^. 

Much fewer studies have been devoted so far to the same questions, namely the temporal 
evolution of two-time quantities, for the case of conserved order parameter dynamics. Re- 
stricting to a system of discrete spins, the rules of Kawasaki dynamics [24] now consist in 
choosing two adjacent opposite spins, and exchange them with a rate depending on the en- 
ergy difference between the initial and final configurations. In the recent past the question of 
the long-time behaviour of the autocorrelation function for conserved dynamics has already 
been addressed (25! 26 , Is^l - I n particular predictions have been given for the values of the 
autocorrelation exponents A and A c governing the decay of the autocorrelation function at 
large temporal separations, respectively in the low temperature phase, and at criticality. Fi- 
nally, a very recent work addresses the question of the response for the Ising-Kawasaki chain 
in the low temperature scaling regime 2^|. The prediction made in this reference states 
that the fluctuation-dissipation plot (that is the relationship between integrated response 
and correlation) of this case is identical to that of the Glauber non-conserved case, which is 
itself known analytically 0,0]. This result is rather surprising because it would imply some 
kind of "super universal" behaviour in the relationship between correlation and response in 
the asymptotic regime. 

The aim of the present work is to revisit and extend these former studies. We investigate 
the behaviour of the two-time correlation and response functions for an Ising spin system, 
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both in one and two dimensions, with Hamiltonian 



(i) 



where (z, are nearest neighbours, evolving under Kawasaki dynamics after a quench of the 
system from a disordered high temperature initial condition to the critical temperature. 

We first give the method used in this paper in order to compute the linear response. 
This method, due to Chatelain j^], and later on clarified in [jy], will be used both for the 
one-dimensional and two-dimensional cases. We then describe the rules of an accelerated 
dynamics for the Ising-Kawasaki chain corresponding to the formal limit T — > 0, which 
is nevertheless faithful, i.e. reproduces exactly the results that would be obtained with the 
usual rules of Kawasaki dynamics for vanishingly small temperatures. We finally present the 
results of extensive numerical computations of the autocorrelation and response functions, 
first for the case of a chain, then for the two-dimensional system at criticality. 

II. MEASURING THE LINEAR RESPONSE WITHOUT APPLYING A FIELD 
IN KAWASAKI DYNAMICS 

In this section we give an analytical expression of the response function to an infinitesimal 
field, for a ferromagnetic system evolving under Kawasaki dynamics, and show how this 
quantity can be measured. We follow the lines of reasoning of refs. jsfll Isij ]. which are 
devoted to the same question for single-spin flip dynamics. 

A. Kawasaki dynamics with heat bath rule 

Hereafter, time t is discrete and counts the number of spin exchange attempts, and not 
the number of Monte Carlo sweeps. In order to define a response function, an external 
perturbing field hi is applied on any site i, and, in presence of the field, the Hamiltonian is 
changed to Ti. — ^2 hi<ii. 

Kawasaki rules consist in updating a pair of two opposite adjacent spins <jj = — (jj (we 
will always omit the S^-a factor in the following), with heat bath probabilities 
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where (3 is the inverse temperature, and the Weiss field hjj takes into account the effect 
of neighbours on the couple of spins which are updated. For a generic 2-spins interaction 
Hamiltonian we have 

kedi\j ledj\i 

where di\j represents the set of neighbours of i, with j excluded. For example, in one 
dimension, with Hamiltonian ((H), we have hj{ +1 = cr^-l — a i+ 2- 



B. Response function 



Following strictly the notation of ref. 



we consider systems made of iV Ising spins, 



where the autocorrelation and the response functions are defined as 



JV N 
i=l i=l 



9±m 

N ^ dhi{a) 



(4) 



with ( • ) representing the average over thermal histories. 

We concentrate on the integrated response function, or susceptibility 



x(t,s)=T / duR(t,u 



(5) 



where the temperature T has been added in the usual definition in order to simplify the 
notation and to have a well defined expression in the T —>■ limit. 

Denoting by I(t) and J(t) the indices of the two spins to be updated at time t, the 
expectation value of the k-th spin at time t is given by 



(cfc(*)) = Tr^) 



°k{t) IJ Wi {t >)j{t>)(p{t')W -I) 



t'=i 



(6) 



where o is the vector of the N spins configuration, the trace is over all the histories a(t' 
with 1 < t' < t, and the transition probability are given by 

exp[/3cr;0^ + K- hj 



IP 



(7) 



2cosh[p(hJf + h i -h j )) 

Note that hjj(a) = hjj(r), because the Weiss field does not depend on the value of spins at 
sites % and j. Since the transition probability only depends on the perturbing fields on 
sites i and j, one has 



dh k 



h=0 



(8) 



where we have defined = tanh(/3/i^). 

Now, if on site k an infinitesimal probing field h k is switched on at time s (i.e. h k (t) = 
h 6(t — s)), all transition probabilities with index k (and only these ones) will depend on the 
perturbing field for times larger than s. Then differentiation of eq. (jUJ) with respect to this 
field yields the integrated response 



Xik{t,s) = T 



d(cri{t)) 



dhi 



h=0 



^i{t)Y[W m ,j {t/) (a{t')\a{t'-l 



t'=i 



[ S I{u),k (^(w) - °Y(u)J(u)( u )) + 5j(u),k (o-fc(w) + (?Y(u)J(u)(u))] 



u=s+l 



(9) 



which can be simply written as a correlation function 



Xik{t,s) = (ai(t) Aa k (t,s)) 



(10) 



where 



Aa k (t,s) = ^ S I{u))k a k (u) - of {u)J{u) {u) + Sj {u ),k cr k (u) + af {u) j {u) (u 



u=s+l 



C. Computing the linear response in a simulation 



Let us note that calculating the linear response in a numerical simulation with no per- 
turbing field using eq. (jlOJ) is as easy as measuring a correlation function. One has to keep 
track of the vector Aa k (t, s) that should be updated for all times between s and t (i.e. when 
the "ghost" field is switched on). At each of these times, one has to: 



compute hjj and a]J = tanh(/?/^ 
update spins and aj according to the heat bath probability (with no external field) 



V{oi = a, Oj 



-a 



2 cosh(pahJJ) 



increment A<jj, A<jj as 



A(jj — > AcTj + (Tj — a 



A 



Co 



A 



°3 + 



W 

ij ! 
W 



(12) 
(13) 
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It is important to note at this point that contributions to the increment Ao~i come either 
when the updated spin flips, or when it does not, keeping its previous value. This is the 
reason why one can speak of a "ghost" field in this method (see above). An illustration of 
this fact is encountered in section 11111 B below. 

The use of eq. ()1L)|) allows one to compute, in a single run, the integrated response (JSJ) 
(zero- field cooled magnetization), as well as the thermoremanent magnetization, 

p(t,s) = T [ duR(t,u), (14) 
Jo 

both for many different values of s. 

III. ACCELERATED DYNAMICS OF THE ISING-KAWASAKI CHAIN IN THE 
ZERO TEMPERATURE LIMIT 

It is well known that Kawasaki dynamics at zero temperature rapidly brings the one- 
dimensional system to a blocked state, where the distance between any couple of domain 
walls is at least In such a situation any spin-exchange move would cost an energy 

AE = 4, and is therefore forbidden. For example TTTLLLII ~~ > TTTITIII would create two 
domain walls. Such a process is called an evaporation 32]. 

However, at an infinitesimal temperature T = e, the evolution may eventually undergo 
such an evaporation with probability 1/(1 + e 4//£ ), hence on a time scale 0(e 4//£ ), which 
diverges for e —>■ 0. After each evaporation process the dynamics proceeds by diffusion, that 
is by moves with energy cost AE = (e.g. |||TIIII — > IIIITIII), until a condensation 
process occurs, corresponding to a move with energy cost AE = —4, which brings the system 
back to a new blocked state (e.g. HITITTT — * HlltttT)- The diffusion and condensation 
processes take place on scales of time 0(1), i.e. much smaller than the typical time between 
two consecutive evaporations, and therefore can be considered as instantaneous if time is 
counted in units of r = e 4,/e , when s — > 0. 

We exploit the strong separation of time scales between evaporation and diffusion/con- 
densation processes, and we simulate the dynamics of the Ising-Kawasaki chain in the limit 
T = e —* with the following accelerated dynamics. Time in this accelerated dynamics is 
counted in units of r = e 4 / £ . 
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A. Rules for the accelerated dynamics 



As long as the system is not in a blocked state, we use the usual T = Kawasaki dynamics 
where only spin exchanges with AE < are accepted. During these lengths of time, since 
the only processes occurring, diffusion and condensation, are actually instantaneous on the 
scale r = e 4//£ , time is not increased at all. 

When the system is in a blocked state (with a number n of domain walls), we choose 
randomly a domain wall, we exchange the two spins on the sides of the domain wall (evap- 
oration), and we increase the time by At = 1/n. 

The choice for At can be understood as follows. In a Monte Carlo simulation at T = e, 
the probability of accepting an evaporation process being approximately equal to e _4/,£ , on 
average, a number 0(e 4/,£ ) of tries will be necessary before a success. In each Monte Carlo 
sweep (MCS) n tries are made, since the only couples of spins which satisfy the requirement 
a", = — CTj are those around a domain wall. Thus the typical number of MCS done before an 
evaporation process takes place is r/n, i.e. At = 1/n in units of r. 

B. Rules for updating Act during the accelerated dynamics 

Being a quantity integrated over time, Acr&(t, s) gets contributions from both fast and 
slow processes. 

Moreover, for T — > 0, eqs. (fP21 H3|) can be simplified to Auj — > A<jj + {2<7j, <jj, 0}, and 
Adj — > A(jj + {2<Tj, <jj, 0}, for evaporation, diffusion and condensation processes, respectively. 

Between two evaporation processes, there is however an additional contribution to take 
into account (see eq. (fTT|) below), which is not apparent in eqs. (fT2l [T3J). This is a direct 
consequence of the fact, mentioned above, that A<Tk(t, s) gets contributions of updated spins, 
even when they do not flip. Consider a blocked configuration, where all the spins are aligned 
with their local fields, <7i = sign(/iJJ) and dj = — sign(Zi^). Since the time spent in each 
blocked state becomes infinite for e — > 0, a non trivial limit for the integrated quantity Act 
is generated. Indeed, working at temperature T = e, and, noting that in a blocked state 
hJI = 2<7j, we have, 




taoh(/ig7e) 



sign(Zi^) 1 - 2e 



*\h%\/e 




) 



(15) 
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hence 



Or - O, 



. w 



2a, : 



10 T 



(16) 



Then the summation over r/n MCS of the last quantity gives a finite limit when e — > and 
t — > oo. In practice this means that, just before leaving a blocked state by an evaporation 
process, all the Act* for spins close to a domain wall must be updated with the following rule 

Act, - Act, + ^. (17) 
n 

Some final comments on the accelerated dynamics described in this section are in order. 

• This dynamics is faithful, i.e. it reproduces exactly the results that would be obtained 
by a standard Monte Carlo simulation at finite temperature T = e, in the limit 
e — > 0. This equivalence will be illustrated below on the example of the mean length 
of domains. It is also very efficient, since it requires much less computational efforts 
than the standard Monte Carlo dynamics. 

• Its definition can be extended to any spatial dimension. 

• Finally, one may wonder how this dynamics compares to the effective dynamics of 
refs. 22, H |- which is only defined in one dimension. The spirit of the later is to 
trace upon all events occurring between the instant of time a spin detaches from a 
domain and that when it reaches the neighbouring domain. The accelerated dynamics 
introduced in this section does not trace upon these events. However doing so -in one 
dimension only- would lead to the dynamics of [3, Em| , and allow a faster computation 
of the average domain length L(t), and with little more work, of the autocorrelation. 
In contrast, tracing upon these events is much more subtle for the computation of the 
response, and would deserve further study. Otherwise stated, it is not clear to us for 



the time being whether the method of 



32, 



34j can be used for the computation of the 



response. The difficulty comes from the fact that Act gets contributions even when 
spins are updated without changing their value. 



S 



IV. RESULTS FOR THE DYNAMICS OF THE ISING-KAWASAKI CHAIN IN 
THE ZERO-TEMPERATURE LIMIT 

In this section we report the results of extensive numerical simulations, using the methods 
described in the previous sections. We are interested in the behaviour of observables in the 
low-temperature scaling regime defined by 1 <C t <C t eq , where the equilibration time t eq is 
related, at inverse temperature (3 = 1/e, to the equilibrium correlation length £ « e 2/3 /2 by 
t cq ~ £ 5 [35]. This regime is naturally attained by the accelerated dynamics. 

A. Mean domain length 



500 r 

450 - 
400 - 




FIG. 1: Mean domain length L(t) against i 1//3 for the Ising-Kawasaki chain obtained by the accel- 
erated dynamics for vanishingly small temperature. Inset: same data on a log-log scale. 



The mean domain length is related to the energy E(t) of the chain at time t by 

Lit) = 2 — . (18) 

W 1 + E(t) K ' 

It is well known that, in the low-temperature scaling regime, L(t) scales as 

m 1/3 

Lit) ~ - , (19) 
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# evap 





10 3 


10 4 


10 5 


10 6 


10 7 


s value 


0+ 


1.1710(5) 


15.998(12) 


312.43(23) 


8000(9) 


235580(740) 


# samp 


895 


1189 


717 


1146 


1179 


51 



TABLE I: For each choice of the number of evaporation processes done before the measurement of 
the autocorrelation and response, the corresponding value of the waiting time s and the number 
of different thermal histories are reported. 



where r = e 4 ^ ~ £ 2 . This scaling is illustrated in ref. |35| . where for increasing values of r 
(i.e. decreasing temperatures) a linear master curve is found when L(t) is plotted against 
(t/r)V3. 

We measured L(t) using the accelerated dynamics of section ILTTl Since time is measured 
in units of r in such a scheme, the curve thus found is the limiting master curve that would 



have been found by the conventional means of |35( in the limit T — > 0. We indeed checked 
that the amplitude of the law (flH|) . L{t) = A (t/r) 1 / 3 , found in [3| was in agreement, though 
less precise, than that found in the present work, A ~ 2.29. Figure Q depicts the master 
curve thus obtained in the T — > limit. The system size is N = 2 12 = 4096, and the number 
of samples are reported in Table |1J (We also checked for the absence of any finite-size effects 
by simulating few iV = 2 16 samples.) 

The inset of figure Q shows the growth of L(t) in a log-log scale. The convergence to the 
slope 1/3 is slow, mainly because of the presence of an offset at initial times. This offset 
corresponds to the first blocked state reached by the system. Remind that the system is 
prepared at time t = in a random configuration. It is then evolved under T = Kawasaki 
dynamics, until it reaches the first blocked state at time t = + in units of r. This first 
blocked state has a mean domain length L(0 + ) = 4.15886(15). 



B. Autocorrelation 



A well-known fact of the kinetics of coarsening with non-conserved dynamics is that, in 
the low-temperature scaling regime, and for large temporal separations (1 <C s ^ t t eq ), 
the two-time autocorrelation C(t, s) decays as |l| 

C(t,s)~(±P-Y\ (20) 



L(s) 
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and, for the particular case where s = 0, as 

C(t,0)~L(i)- A , (21) 



defining the autocorrelation exponent A 

The case of conserved dynamics is more complex. For quenches to temperatures below 
T c , Yeung et al. [28] find bounds on the autocorrelation exponent which depend on the value 
of the smallest time s. For s — 0, A > d/2, where d is the dimensionality of space, while for 
values of s in the scaling regime, A > d/2 + 2 for d > 2, and A > 3/2 for d — 1. 

Given the prediction of ref. [25J that, in the low-temperature phase, C(t, 0) ~ L(t)~ x , 
with A = d, Yeung et al. conclude that, for d = 1 (and more generally for low dimensions), 
the behaviour (|21|) holds for s small, while for s in the scaling regime the behaviour ()20|) 
should be replaced by 

L(t) - x 



c(t - s) ~UMJ (22) 

with a different exponent A' > A. 

The prediction above therefore implies that for the Ising-Kawasaki chain at vanishingly 
small temperature the curve of the autocorrelation for two values of s, one being taken in 
the short-time regime, the other one in the scaling regime, should cross at some later time 
t*. This simple observation was not noticed by the authors of The question therefore 
arises of what is the behaviour of the autocorrelation C(t, s) for later times t^> t* . 

Hereafter we suggest the following scenario, based on the reasonable hypothesis that the 
two curves mentioned above do not cross. Define more precisely t*(s) by 

i.e. t*(s) ~ s A 7( A '- A ) # Then, there are two scaling regimes: 

• the intermediate scaling regime of Yeung et al., for f<f, where holds, 

• the ultimate scaling regime, for t ^> £*, where 

C(t,s)~L(t)~\ (24) 

Therefore the following scaling law should hold 

C(t,s)=L(t)-'g(^pi, (25) 



L( S )* 
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with g(x — > 0) ~ x 1 in the intermediate regime, while in the ultimate regime g(x — > oo) 
should converge to a constant. 




10 100 

L(t) 



FIG. 2: Two-time autocorrelation function of the Ising-Kawasaki chain at vanishingly small tem- 
perature. Full line has slope —A = —1, while dotted ones have slope —A' = —2.5. 



Outcomes from our simulations are compatible with these predictions. Figure El shows 
the autocorrelation function as a function of L(t), for the values of s given in Table UJ The 
presence of two different scaling regimes is evident, although the intermediate regime is not 
very clean on small s data, and the ultimate regime is not reached for large s data. Figure El 
is a plot of the scaling function g(x) defined above. We find A' ~ 2.5, with A = 1. 

Let us note that before attaining the intermediate scaling regime there is yet another 
temporal regime, clearly visible on figures El and El which takes place at very short times. 
One can show that in the scaling variable of figure El this very short time regime becomes 
smaller when s increases, and therefore irrelevant in the s — > oo limit. Indeed one finds 
numerically -see figure 0]- that in the very short time regime the autocorrelation function 
can be written as follows 

1-CM«^ • (26) 
Eq. (J2T?|) implies that the timescale of decorrelation processes happening in the very short 
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5 = 235580(740) ' — ' — ' 




s = 8000(9) 




5 = 312.43(23) 




5= 15.998(12) ■• a i 




5=1.1710(5) ------ 




v = ()' 


1=1 V =2.5 ^1 








0.01 0.1 


1 10 



L(t) X '~ X / L(s) X ' 



FIG. 3: Scaling function for the two-time autocorrelation of the Ising-Kawasaki chain at vanishingly 
small temperature (see text). 



time regime is s 2 / 3 . This timescale is much smaller than the timescale of the intermediate 
regime, which is s. 

In preparation of section |VJ let us mention that, at criticality, dynamical scaling predicts 
that 

C(t,0) ~ L(t)~ x % (27) 

defining the critical autocorrelation exponent A c . We have also, for s in the scaling 
regime, 

C(M)~L(,)*(|)^ (28) 

where (3 and v are the usual static critical exponents (/3 = 1/4, v = 1, in 2D). This form 
holds for both non-conserved and conserved dynamics, with the possibility that, for the 
latter, the exponent appearing in (|28j) be not the same as in (|27|) . as discussed in section IV! 
For conserved dynamics at criticality, refs. (25), 26, predict A c = d, while statements 
made in 28] on the long-time behaviour of the autocorrelation are less precise. 

A last comment is in order. At criticality, there are no well-defined growing domains. 
The interpretation of the growing length in eqs. (}2Tj) and (|2*Hj) is the typical size over which 
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the system looks critical. The present situation of a system evolving after a quench from 
high temperature to T = T c is nevertheless usually referred to as "critical coarsening" . By 
convenience we shall still call the length L(t) the mean domain size. 

Note that for the Ising chain the magnetization exponent (3 = 0, hence there is no 
distinction to be made between the two behaviours (|2L))) and (|28|). 



C. Response and fluctuation-dissipation plot 



Let us recall that for the zero-temperature non-conserved (Glauber) case, the two-time 
correlation and response functions are known analytically ^, As a consequence, the 
fluctuation-dissipation ratio X(t,s), defined by [5( 

m «\ - x ^ s) dc ^ s) 

can be obtained in closed form. In the scaling regime, it reads 

*(*,.) + 

yielding, for large temporal separations, the limiting ratio 

= lim lim X(t,s) = -. 

s^oo t^oo 2 
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In the scaling regime, the integrated response functions p(t, s) and s) only depend on 
t/s, or equivalently on C, 

p{C) = ±= - X (C) = ^ arctan (±= tan ^ . (29) 

Alternatively the limiting ratio can be extracted from p(C), as 

.^=Hm^f>. (30) 

Note that, in contrast to the generic two-dimensional case considered in section for the 
Ising chain, the integrated response function does not bear any dependence in s because the 
exponent (3 = 0. 



?3 



s = 8000(9) 

5 = 312.43(23) 

5= 15.998(12) 

5= 1.17053(58) 
y = + 

Glauber 




FIG. 5: Integrated response against autocorrelation for the Ising-Kawasaki chain at vanishingly 
small temperature. 



We now turn to the Ising-Kawasaki chain. Figure El depicts a plot of the integrated 
response x{t, s) against the correlation C(t, s), in the limit of vanishingly small temperature, 
for the values of the waiting time s given in Table HJ Most of the plot corresponds to the 
intermediate regime (t t*(s)), since the ultimate regime is reached only for the two 
smallest s values and very small correlations C(t, s) < 0.03 (see figure^- Restricting to the 
first regime, we observe a slow convergence of the data to a limiting curve when s increases. 
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However this limiting curve lies above the theoretical Glauber curve (|29j) . in contradiction 
with the prediction made in ref. 2^|, stating that the fluctuation-dissipation plot for the 
Ising-Kawasaki chain is identical to that of the Ising-Glauber chain, for vanishingly small 
temperature. 

We are thus led to critically review ref. j^]. A similar plot of the integrated response 
against correlation is presented in this reference, for three sets of values of T and s: (T = 
0.48, s = 3 • 104), (T = 0.70, s = 2 • 103), and (T = 0.70, s = 104). For these values of (T, s) 
the data follow rather closely, at least in a range of values of the correlation, the theoretical 

n 

Glauber curve ()29j) . The authors of ref. [29] therefore induce that the data obtained for 
vanishingly small temperatures and increasing values of s should eventually converge to 
eq. J2D. 

Noting that the values of (T, s) mentioned above correspond, in units of e 4 ^ 3 , to s — 7.2, 
s = 6.6, and s = 33 respectively, we see on figure El that, indeed, for this range of values of 
s, the data points fall not too far away from the Glauber curve. However, since this holds 
neither for smaller nor for larger values of s (in units of e 4 ^), we conclude that the apparent 
identity between Kawasaki and Glauber curves observed in [^J is just a coincidence due to 
the range of values considered in this reference. 

The existence, in the intermediate regime, of a non trivial scaling limit for the response 
implies, as a corollary, a non trivial value of X^, when the ratio x = t/s —>■ oo. A precise 
numerical value of this quantity is however difficult to obtain from ZFC data. At present 
we can not exclude the possibility that is the same for both Glauber and Kawasaki 
dynamics. 

Consider now the ultimate regime (t ^> t*(s)). This regime corresponds to C < C*(s) 
where C*(s) = C(t*(s), s) ~ S -^7(A'-A) ) with AA'/(A'-A) ~ 1.66. This is reminiscent of the 
situation encountered at criticality for the two-dimensional Ising-Glauber system [10] . In this 
case the fluctuation-dissipation theorem holds except in a region C < C*(s) ~ s~ ac , where 
a c = 2f3/uz c ~ 0.115, which vanishes for increasing values of s. However this mechanism 
does not prevent the occurrence of a non trivial limiting ratio X^ (nj. (See section El) By 
analogy, we expect a non trivial value of X^ in the ultimate regime, a priori different from 
that obtained in the intermediate regime. Note however that since the exponent AA'/(A' — A) 
is larger than 1, C*(s) is decreasing very fast, and, as a consequence, the regime where X^ 
could be measured is hardly reachable in a numerical simulation. 
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V. RESULTS FOR THE CRITICAL DYNAMICS OF THE 2D ISING-KAWASAKI 
MODEL 

The aim of this section is to investigate the critical coarsening of a two-dimensional system 
of spins evolving under Kawasaki dynamics from a random initial condition. We will, as in 
the previous section, and inspired by the results found there, examine in turn the behaviour 
of the two-time correlation function, then of the two-time response function. 



A. Mean domain length 



The numerical study of the long-time behaviour of C(t , s = 0(1)) for the critical dynamics 



of the 2D Ising-Kawasaki model is well established 
holds, with a decay exponent A c = 2, confirming the prediction o 
itself is observed to grow as L{t) ~ t 1 ' Zc , with z c = 4 — 77 = 15/4 



261 ] . Scali ng i s observed, i.e. eq. (|2~7jl 



3 



25| . The mean domain size 
A confirmation of these 



results is provided by figure El where the mean size of domains is extracted by the excess 
energy with respect to the equilibrium energy at T c , E eq = l/v2, i.e. L(t) = (E(t) — E ec Y^. 
This method has been alread y u sed for low temperature coarsening, either for conserved (37 1, 



or non conserved dynamics 



361 ] . It is however the first time that it is used at criticality, 



where in general the growing length scale L(t) is obtained from the position of the first 
zero of the equal-time correlation function. Our reason to do so lies in the fact that for 
critical coarsening dynamical scaling holds, as it does for low temperature coarsening, and 
that therefore there exists only one single growing length scale in the system. 

In order to compare the qualitative behaviour of both conserved and non-conserved dy- 
namics at T c , we take a series of snapshots at instants of times where similar domain sizes 
were reached in both dynamics (see figure EJ). Since for non-conserved dynamics z c « 2.17, 
in 2D, conserved dynamics is much slower. It is interesting to note the overall similarity of 
the snapshots at corresponding instants of time. 



B. Autocorrelation 



The study of the behaviour of C (t, s) when the waiting time s is deep in the scaling 
regime is largely unexplored. We first measured the autocorrelation as a function of L(t) 
for different values of s (see fig. (HJ). As in the one-dimensional case, we observe a seemingly 
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1 10 

L(t) 

FIG. 6: 2D Ising-Kawasaki model at criticality: C(t, 0) and L(t). 




t=20 t=93 (=41122 t=551247 



FIG. 7: Snapshots of configurations after a quench from T = oo to T = T c in a 500 2 Ising spin 
system. The four snapshots on the left correspond to Glauber dynamics, while the four on the 
right to Kawasaki dynamics. Times have been chosen in order to have similar values for L(t). 
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different decay exponent as soon as s is large enough, i.e. in this regime eq. (J2EJ) should be 
read with a different exponent A^.. 




FIG. 8: 2D Ising-Kawasaki model at criticality: C(t,s) against L(t) for different values of s. 



In order to assess this point, we define, by analogy with the one-dimensional case (see 
eq. (|25|)). the scaling function g c as 

C(t, i ) = £(t)^,( I ^ ), (31) 

with g c {x — > 0) ~ x^ 1 in the intermediate regime, while in the ultimate regime g c {x — > oo) 
should converge to a constant. Figure |H1 depicts the scaling function obtained using A c = 2 
and = 3.5. In view of this figure it is reasonable to conclude again in favour of the 
existence of these two scaling regimes, defined by the relative magnitude of t with respect 
to the crossover timescale t*(s) ~ s( A c -2 0/ I/ )/( A c-*c). Otherwise stated, t <C t*(s) in the 
intermediate regime, while t ^> t*(s) in the ultimate regime. We have, with the values of 
the exponents given above, and 2(3/ u = 1/2, t*(s) ~ s 2 . 



C. Response and fluctuation-dissipation plot 

We now turn to the response. Following Q], we choose to compute the thermoremanent 
magnetization (fT^j) (TRM). Due to the fast increase of the crossover timescale, t*(s) ~ s 2 , 
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FIG. 9: 2D Ising-Kawasaki model at criticality: scaling function for the two-time autocorrelation. 

and to the extreme difficulty to have precise data for the response at very long times, the 
results presented below only concern the intermediate regime. These results are contained 
in figures CHI d andd 

In the intermediate regime we assume for the TRM a scaling form similar to that of the 
correlation function, equation ()28|) . that is 



This scaling form is well verified by our numerical data, as illustrated by figure ITUl which de- 
picts plots of the rescaled correlation and response functions, both for Glauber and Kawasaki 
dynamics. For the latter, it is interesting to note that, when s increases, the master curve 
is attained from above for the correlation, and from below for the response, indicating that 
asymptotically the two scaling functions should have the same algebraic decay with exponent 
A'. Plots of the scaling functions of autocorrelation and response for the two-dimensional 

n 

Ising-Glauber model at criticality first appeared in [10j. The value of the exponent measured 
in the present work is slightly smaller than that found in jiol ]. 

Another representation of the same data is given in figure ITU which depicts the parametric 
plot of the rescaled response against the rescaled correlation, both for Glauber and Kawasaki 




(32) 
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FIG. 10: 2D Ising model at criticality: scaling functions of correlation and response, (left) Glauber 
dynamics. The full lines have slope A c = —1.54. (right) Kawasaki dynamics. The dotted lines have 
slope A^, = —3.5. 

dynamics. This figure shows that the two dynamics lead to the same phenomenology. Indeed, 
define a crossover scale for the autocorrelation by C*(s) = C(2s, s) jlOj, corresponding, for s 
large enough, to a value of the abscissa on figure ITT1 approximately equal to 0.47 for Glauber 
dynamics, and to 0.27 for Kawasaki dynamics. Then, on the right part of the plots with 
respect to these values, the slope of the two curves is equal to one, in agreement with the 
fluctuation-dissipation theorem. The left part of the plots corresponds to the scaling regime, 
with a crossover towards a non trivial slope at the origin, equal to the limiting violation ratio 



In order to extract the numerical estimate of these slopes from these data, we plot them 
in two different fashions, as shown on figure IT2l Since p(t, s) ~ C(t, s) (see (fHo"|) ). we plot 
the ratio p/C, first versus L(s)/L(t) (figure a), then versus the ratio of times s/t (figure b). 
Error bars are not shown in the latter, in order to improve readability. For Glauber dynamics 
seems to be attained linearly in s/t , which leads, by taking the intercept with the y- 
axis, to theprediction X^ lauber ss 0.33, in agreement with the estimates given in recent 
works Note that for the spherical model the correction to is exactly in s/t [10]. 

Finally, though it is difficult to be conclusive on the sole basis of figure H21 the latter gives 
more evidence in favour of different values for the limiting ratios corresponding to the two 



X, 



oo • 



dynamics, with X ( 



Kawasaki 




oo 
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C(t,s)L(s) 



FIG. 11: 2D Ising model at criticality: parametric plots of rescaled response versus rescaled corre- 
lation for Glauber and Kawasaki dynamics, using the data of figure ITU1 The full line at the origin 
has slope 0.33. The other one has slope 1, and is meant as a guide to the eye. 




FIG. 12: 2D Ising model at criticality: (left) p/C versus L(s)/L(t); (right) p/C versus s/t. 



VI. CONCLUSION 



Let us summarize the most salient points of this study. 

We extend the method of (30, 31] to the case of Kawasaki dynamics. We introduce a new 
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method for the investigation of the low-temperature scaling regime of the Ising-Kawasaki 
chain. We define rules for an accelerated dynamics, which are both faithful and efficient. 
This method can be extended to higher dimensions. We find new results concerning the 
behaviour at large time of the autocorrelation function for the critical Kawasaki dynamics 
in both dimensionalities, and demonstrate the existence of an ultimate regime, which was 
overlooked in previous studies. We believe this regime to be also present in the behaviour of 
the response, though hardly accessible with present computer capabilities. As a corollary, 
we expect the existence of a different value of the limiting ratio in the ultimate regime, 
which would require time scales which are presently unreachable. 

In the course of this study we were led to question the validity of the results of ref. 
concerning the fluctuation-dissipation plot for the Ising-Kawasaki chain. The evidence, 
claimed in 2^1 , for the identity of the fluctuation-dissipation plots for Glauber and Kawasaki 
dynamics is seemingly coincidental, and lies in the range of values used in this reference. 

For the two-dimensional system at criticality, the fluctuation-dissipation plots obtained 
from the two dynamics are phenomenologically similar, but not quantitatively identical. In 
particular, the limiting violation ratios are found to be different. It is harder though to 
be conclusive in the two-dimensional case than in the one-dimensional one. In the former, 
the scaling region defined by C(t,s) < C*(s) ~ s -o.ns j s sma ii ; while in the latter, since 
(3 — 0, it identifies to the whole range < C(t, s) < 1. 

In the intermediate scaling regime numerical data are compatible with the simple form 
\' c = d + 3/2 ior the autocorrelation exponent. It would be interesting to assess the validity 
of this hypothesis. A theoretical explanation of the existence of the ultimate regime, both 
for the autocorrelation and for the response, is needed. Finally, a natural extension of 
the present work is the study of the behaviour of the autocorrelation and response in the 
low-temperature phase (T < T c ) of the two-dimensional kinetic Ising model with Kawasaki 
dynamics. 
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